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THE PROBLEM OF INTERNAL AND EDGE CRACKS 
IN AN ORTHOTROPIC STRIP* 

by 

F. Del ale and F. Erdogan 
Lehigh University 


ABSTRACT 

The plane elastostatic problem of internal and edge cracks in an 
infinite orthotropic strip is considered. The problems for the ma- 
terial types I and II are formulated in terms of singular integral 
equations. For the syirmetric case the stress intensity factors are 
calculated and are compared with the isotropic results. The results 
show that because of the dependence of the Fredholm kernels on the 
elastic constants, unlike the crack problem for an infinite plane, 
in the strip the stress intensity factors are dependent on the elastic 
constants and are generally different than the corresponding isotropic 
results. 


*This work was supported by NASA-Langley under the Grant NGR-39-007-011 
and by NSF under the Grant ENG 73-045053 AOl. 


1. 


INTRODUCTION 


In plane elastostatic problems for an infinite orthotropic medium 
containing a line crack [1-3] or a series of collinear cracks [4] it 
was shown that the stress intensity factor is identical to that found 
for isotropic materials. However, if the geometry of the medium is 
bounded, it is expected that in orthotropic solids the material con- 
stants would influence the stress intensity factors. A bounded speci- 
men geometry which is sufficiently simple for the purpose of analysis 
and at the same time is of sufficient practical interest is that of a 
long strip containing internal or edge cracks. The main objective of 
the present paper is by considering this problem to give some idea 
about the degree of influence of the material orthotropy on the stress 

il 

intensity factors. The equivalent isotropic case is one of the more 
widely studied crack problems in technical literature (see, for ex- 
ample, [5-13]). 


2. FORMULATION OF THE PROBLEM 


Consider the plane problem for the orthotropic strip shown in 
Figure 1. Referring to, for example [14], the equilibrium equations 
for an orthotropic plane may be expressed as 
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where u, v are the x,y- components of the displacement vector and 
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for generalized plane stress, and 
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for plane strain. Here, v^, 6.^ are the engineering elastic con- 
stants, (i,j) = (1 ,2,3), the indexes (1,2,3) corresponding to the direc- 
tions (x,y,z) , and the matrix (b. .) is given by 

IJ 


(b^j.) = B = A""* , A = (a^j.) , (i,j) = {l,2,3) , 

a.. - l/Ej. i a^ . = -V../E.. = a.. , (ij*j) 
n ■ ii ’ ij ij'n ji * ' 


(4) 


The solution of the problem shown in Figure 1 may be obtained by using 
the standard superposition technique. Thus, for the purpose of evaluat- 
ing the stress intensity factors and obtaining information relevant to 
the fracture of the solid, it is sufficient to consider the problem in 
which statically self-equilibrating crack surface tractions are the 
only external loads. 


To solve the differential equations (t)let 

u(x,y) = |- J^[fl(a,x)cosay + g^(a,y)s1nax]da , 

v(x,y) = |^[f 2 (a,x)sinay + g 2 (a,y)cosax]dot . 

Substituting from (5) into (1) the functions f^ and g. , (1=1,2) are 
obtained as follows: 


f](atX) = I Aj(a)e ^ 


, f2(a>x) = I CjAj(a)e ^ 




( 6 ) 


s.ay/Ps <♦ s^ay/Ps 

9-\ia»y) ~ f Bj(a)e ^ , 92(a,y) = J djBj(a)e , 

where are the roots of the following characteristic equation: 

s“ + B^s^ +02 =0 , S3 = -s^ , = -S2 . (7) 

and the constants 0^, 0g, Cy and d^, (j«li..,4) are given by 

0^ = (0^-0^02-l)/Pi , = &2/^l ’ 

c^ “ ”C3 ® (1 0iS|)/03Si f C2 “ “C^ “ (1 “ ^^821^352 » ( 8 ) 

d-j “ -dj = “ BiB5)/32Si0g > d2 “ 

Assuming that x and y are axes of symmetry for loading as well as 
geometry, the unknown functions A. (a) and B.(a), (j=l,..i4) are de- 

J J 


termined from the following conditions: 

u(x,y)->0 , v(x,y)*H) for y+« , (9) 

0 , cfj^y(h,y) - 0 , 0<y<“ , (10) 

u(0,y) = 0 , o^y(0,y) = 0 , 0<y<~ , (11) 

a^(x,0) = 0 0<x<h , (12] 

Oyy(x,+0) = -p(x) , a<lxl<b , 


v(x,0) = 0 , 0<|xj<a , b<|xl<h , 


(13 




where the crack surface traction p(x) is a known function. The seven 
horrrogeneous conditions (9-12) may be used to eliminate seven of the un- 
knowns and the mixed boundary conditions (13) would give a system of 
dual integral equations to determine the eighth. In this paper, de- 
fining a new unknown function 

(}>(x) = v(XiO) , a<lxl<b , (14) 

the problem is reduced to a singular integral equation by using the 
first equation of (13). From the second equation of (13) it is seen 
that 

(j)(x,0) = 0 , 0<lx|<a , b<lx|<h , 

fb 

<|.(x)dx « 0 . (15) 

Examining the roots of (7)- it may be observed that 

(i) for 0^<O, 3g = 3|-4$g>0 there are four real roots, s-j, S2 
Sj- -s.| , and s^^ -Sg (s^>0, S2>0); in this case the corresponding 
material will be classified as type I; 

(ii) for 3g<0 the roots are complex, s^j = + iw^® -Sj, 

Sg = -iti)2 = -s^ (w^>0, U2>0) and the related material will be classi- 
fied as type II; and 

(iii) for $^>0, 3g><7the roots would be pure imaginary 

S-| — lUg ® ^2 ~ ^^4” '"^4’ 

In practice generally 0^ is negative. Therefore the problems of 
interest are those relating to materials type I and II only. This 
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classification seems to be necessary In order to pursue the formulation 
of the problem beyond equations (6) without Introducing unnecessarily 
complicated complex algebra. Also, 1n this paper only the case of 
generalized plane stress will be considered. For plane strain It Is 
sufficient to replace the quantities • V 

VyxEj(/(l-VxyVyx) by b^^, b 22 , and b^g* respectively (see equation 4). 

Because of symmetry considering only one quarter of the medium 
shown In Figure 1, and using the standard stress-displacement relations 
for plane stress, after somewhat lengthy but routine analysis, for ma- 
terial type I (i.e., for real and Sg) the problem may be reduced to 
the following Integral equation: 


I Ct^+:^+k^{x,t) - ki(x,-t)]4»(t)dt 

7t(1-V V ) 

under the additional condition (15). Here the kernel is given by 

-(b-t)aB5/s^ 


(16) 




mi 4 

•’(h“t)aB5/s2 

+ K2(x,a)e ]da . 


(17) 


The expressions for K.j, Kg, and m^^ are given in Appendix A. 

For material type II the roots of the characteristic equation (7) 
are complex. Defining 

2= -S3 , S2 = oj^ -1(1)2= -S4 ’ 


S| = w-j+iu, 


(18) 


and, assuming that (i)<|>0, in this case the integral equation becomes 




* ' ' 2? r ^ » 

y*^i4 

again subject to condition (15). The kernel is given by 

-u)^a(h-t) 


(19) 


“ Iq K2(x,t,ct)e 


da , 


( 20 ) 


where the function and the related constants r^ are defined in the 
Appendix B. 


One may note that in the special case of single internal crack* 
(i.e., for a=0, b<h) the integral equations (16) and (19) may be 
written as 



(x,t)]({)(t)dt 


p(j^) » "b<x<b , 


(i=l,2) 


*^1 




xy'^yx 






V 

xy yx 


( 21 ) 


where i=l and 2 correspond to materials type I and II, respectively. 
In this case the single-valuedness condition (15) becomes 


|-b 

i'(|)(x)dx = 0 . 

4 


( 22 ) 
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3. STRESS INTENSITY FACTORS 

The standard definition of the stress intensity factors at the 
crack tips a and b is 


k(a) » lim *^Z(a-x) a^^(x.O) , 


k(b) » lim »/2(x-b) a (x,0) . 

x-^b yy 


(23) 


To determine the asymptotic behavior of the cleavage stress a around 

*/ «/ 

the crack tips, first it may be observed that the index of the singular 
integral equations (16) and (19) is +1. Therefore, the solution is of 
the following form [15]: 

-h 


<|)Vt) » f(t)[(t-a)(b-t)] 


(24) 


Next, it should be pointed out that the left hand side of (16) and (19) 
gives cfyy(x»0) for x outside the interval (a,b) as well as within. 

Thus, making a change in variable 




(25) 


for example, (16) may be expressed as 




where 


(26) 


q(s) = Oyy(x,0) , ij^(r) = <l)(t) = F(r)(l-r^) 


(27) 
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and k(s,r) and F(r) are bounded functions. The objective is then to 
determine the asymptotic behavior of q(s) around s=+l, (|s|>l) in terms 
of the unknown function F(r) which U. obtained by solving the integral 
equation (26) in -l<s<l where q(s)=-p(x) is known. To do this let us 
assume that the bounded function F(r) can be represented by the follow- 
ing infinite series in Tchebyshev polynomials: 

F(r) = f T„(r) . (28) 

0 

Substituting from (28) into (26) one obtains 


|A„CG„(s).Hn(s)] = ^ 
S (s) - i f’ 

" ^ ll (r-s)/!^ * 


Hn(s) 


, rl k(s,r)T (r) 
i 0 — dr 


(29) 

(30) 


Here, H (s) is a bounded function. To determine G„ one may use the 
’ n n 

expression 

1 [' , ( z-^" (3,) 

^ il [r-z)/F^ 

where z is the complex variable in the plane cut along (-1 ,1) . Ob- 
serving that on the real line (z^-l)^ is an odd function, from (30) 
and (31) it follows that 

6 (s) . - L s-S9n(s) . ./FT 1.. , (32) 

" sgn(s).^^ 

As s-»+l (32) yields 


- 9 - 


(33) 


s-H-1, (s>l): G|^(s) = - R,(s) , 

s»-l, (s<-l): G^(s) • R 2 (s) . 


where the functions and Rg are bounded at s=+l. 

Now, observing that » finite, 1^,(1) » 1, Tj,(-1) » (-1)”, from 

(29) and (33) the asymptotic behavior of q(s) is found to be 

(S>1): .^+R3(s) . 

S+-1, (s<-l): ^ » -^+R,(s) . (34) 


where again the functions R^ and R^ are bounded at s»tV. Going back 
to the original quantities by using (25) and (27), (34) becomes 


x->b. (x>b); + Rc(x) 


x-^a, (x<a): 


Ml 5 

(x,0) 


M 


1 


• s# • ■.<-) 


(35) 


where the functions Rg and Rg are also bounded at x=b and x=a. Thus, 
from (23) and (35) the stress intensity factors are found to be 

k(b) = - M^F(l)/t , k(a) = M^F(-l).^ . (36) 

In the case of fully imbedded cracks the integral equations (16), 
(19) or (21) can always be reduced to the normalized form (26) and can 
be solved by using the technique described in [16]. 
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4. EDGE CRACKS 


In equations (16) and (19) the kernels k-[(x,t)-k^(x»-t) and 
k 2 (x,t)-k 2 {x,~t) are bounded provided b<h (see Figure 1), For b=h, 
that is in the case of edge cracks, the integral equations are still 
valid but these kernels do not renrain bounded as x and t go to the 
end point b»h ana, consequently, the singular behavior of the solu- 
tion at X“b=h wuy no longer be described by (24). Expressing the 
kernels in (16) and (19) as 

k^(x.t) = k.^(x,t) + k^g(x,t) , (i»1.2) 

where is bounded in the closed interval [a,h]» the unbounded parts 
k^g, (i=l,2) may be obtained from (17) and (20) by examining the asymp 
totic behavior of the integrals for large values of a. Thus, after 
some routine analysis we find 







( 38 ) 


r2QW^{2h«t-x) 
wi(2h-t-x)^ + w^t^h-t-x)^ 

where the constants and are given in the appendices. Thus, for 
example, the integral equation (16) may be expressed as 

j " k^(x/t)]«j>(t)dt 

= - J p(x) , a<x<h ♦ (39) 

In (39) for the purpose of asymptotic analysis transferring the 
terms involving the bounded kernels to the righthand side one may 
write 

fh , 

where is a bounded function in [a,h]. Letting now 
f,(t) 

^(t)=_-J_ ^ , 0<Re(a,8)<l , - (41) 

(h-t)“(t-a)P 

where f^ is H-continuous in [a,h], and following the procedure outlined 
in [15], the characteristic equations for u and 8 are found to be 

cotTTg = 0 , 3*1/2. (42) 

-cosTia + [m^g(3s/s?)“(si/35) + 111 ^ 7 ( 51 / 35 ) 

14 15 

+ nilQ(s2/3s) + m^g(35/s2)“(s2/85)] = 0 . (43) 
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Similarly* for material type II described by (19), assuming the 
solution again as given by (41), the characteristic equations become 


cotirp * 0 , 3 = 1/2 


( 44 ) 


COSTTa + 


"’14"19 


’^29 ^ -“ISfTSr' 


% 


“ uT+"wf sin(2atan jj^) 


(45) 


At the imbedded crack tip x=a it is seen that the singularity has the 
expected 1/2 power. On the other hand, as in the isotropic case, (43) 
and (45) have no root for which 0<Re(a)<l, meaning that at x=h there is 
no power singularity. One may also proceed and investigate the possi- 
bility of a logarithmic singularity for the solution. Thus, letting 
a=0 in (41) and defining the sectionally holomorphic function 




dt 


(46) 


we find [15] 


fn(a)e' 


,ffi3 


fi(h) 


1 iiX"/ 

• 'sihVe '• 7 ' "sS ' ~T — :b * 

1 simr3 (2_a)3 TT(h-a)^ ^ 

f ' 

1 A(t\ fi(a)cotiT3 f,(h) 


J L i-l^h-xT " •^i(2h-x) , ... 


(47) 


- 13 - 


where ?2 ^3 bounded at x=h and have at most a singularity of 

lower order than 3 at x^a^ Substituting from (47) Into (40), multi- 
plying through by (x-a)^ and letting x-»-a, it Is found that cot7rp=0, 
giving again 3=1/2. On the other hand. In the neighborhood of the 
end point x=h one obtains 


[ 1 - 


”’l4"’l5^5 


(s^m^g+s^m^^+S2m^g+S2m^g)]log(h-x) + P^(x) = P^(x) (48) 


where P^ contains all the bounded terms around x=h on the left hand 
side of (40). Similarly, for the material type II one finds 


(49) 

\i 

It turns out that, as in the case of isotropic materials [10}, the co- 
efficient of the logarithmic term In (48) and (49) is identically zero, 
meaning that the solution may not have logarithmic singularity at x=h. 
In the edge crack problem th^ integral equatici? 19) and the similar 
equation for the material type II are solved by defining 


^(x) 


Kxl = Hil 

/s+1 


(50) 


and by using the numerical technique described In [10]. In this case 
the stress Intensity factor at x=a becomes 

k(a) = M,F(-1)*^ , A = (h-a)/2 . 


(51) 



5. RESULTS AND DISCUSSION 


As an example the following two orthotropic materials will be 
considered: 

Type I: E^ = SxlO^ psi (55.16x10® N/m2) , 

Ey =^2 h.75x10® psi (170.65 x109 N/ro2) , 

G = 0.7x10® psi (4.83 xl09 N/m^) , 

AJ 

« 0.036 , 

Type II; E^ »^ 3.1 x 10® psi (21.37 x1q 9 N/m^) , 

E = 9.7x10® psi (66.88x 109 N/m^) , 

W 

G^„ = 2.6x 106 psi (17.93x109 N/m^) » 

xy 

Vxy = 0.2 . 

Tables 1-3 show some of the calculated results for the stress inten- 
sity factors. In all the calculations it was assumed that the crack 
surface traction was constant, i.e., 

ffyy(x.O) = -p(x) = -Pq (52 

which corresponds to uniform tension of the strip away from the crack 
region. Table 1 shows the results for an internal crack of length 2b 
(see Figure 1, a=0) which was found by solving (21). The stress in- 
tensity factors used in the tables are defined by (23) and are calcu- 
lated by using (36) for internal and (51) for edge cracks. The stress 
intensity factors for symmetrically located two collinear cracks 
(Figure 1) are given by Table t. Table 3 gives the results for sym- 
inetric edge cracks. 



The tables also contain the stress Intensity factors for the 
Isotropic strip which are included for comparison. A close examina- 
tion of the integral equations (16), (19), or (21) would indicate that 
in orthotropic materials since the Fredholm kernel k-j or k2 Is heavily 
dependent on the material constants, the solution must also depend on 
the constants. On the other hand, in isotropic materials even though 
the structure of the integral equation is identical to that of (16) or 
(19) (see, for example, [10]), the kernel of the integral equation is 
independent of the elastic constants and the constants appear in the 
equation as a multiplying factor (in the form of (1 +k)/4u) only. The 
stress intensity factors given in the tables indicate that the results 
for the orthotropic strip are indeed different than the isotropic re- 
sults. The tables also show that for approximately the same modulus 
ratio Ey/E (in the example approximately 3/1), depending on the re- 
maining constants, the materials may not only be of different type 
(I or II), also the stress intensity factors may be greater (in this 
case, in material type II) or smaller (in material type I) than the 
isotropic values. In orthotropic materials there are three independent 
material parameters, namely, G /E^, E /E , and v . Therefore, it 
does not seem to be feasible to make a systematic study and demonstrate 
the effect of the material orthotropy on the stress intensity factors. 
However, it appears that there exists a difference between isotropic 
and orthotropic results and in highly orthotropic materials it may 
be significant. 

In solving this problem, the numerical analysis produced a some- 
what unexpected result. First, it should be pointed out that the 



results given In the tables are accurate to roughly three significant 
digits, the remaining digits may not be accurate. On the other hand, 
after rotating the roatevial 90 degrees (i.e., taking the strip parallel 
to the less stiff axis and the crack along the stiffer axis) and fully 
expecting to obtain a different set of results, the print out for the 
stress Intensity factors came out to be Identical - In all eight digits* 
to the original values obtained for the strip which was parallel to the 
stiff axis. Furthermore, the ratio of the function F(r) defined by (27) 
at all points In -l<r<l for the 0 and 90 degree orientations was found 
to be constant, indicating that the crack surface displacements for the 
two cases are related by (see (14), (27), and (36)) 

Vq(x,0)Mo ® (53) 

where the constant M Is defined by (21). This simply shows that the 
kernels k^ and which appear in the Integral equations (16), (19), 
and (21) remain invariant under a 90 degree rotation for a given ortho- 
tropic strip. 


Table 1. The stress intensity factor k(b)/Po‘^ for an Internal crack 
of length 2b in isotropic and orthotropic strips. 


b/h 

Isotropic 

(Jrtlio tropic I 

Tyjf>e I. 

Type II 

-►0 


. 0 

+1.0 

0.1 

1.0060 

1.0044 

1.0064 

0.2 

1.0246 

1.0182 

1.0261 

0.3 

1.0578 

1.0428 

1.0611 

0.4 

1.1094 

1.0811 

1.1155 

0.5 

1.1869 

1.1387 

1.1966 

0.6 

1.3033 

1.2264 

1.3183 

0.7 

1 .4888 

1.3674 

1.5099 

0.8 

1.8160 

1.6241 

1 .8471 

0.9 

2.5809 

2.2487 

2.6278 



Table 2. The stress intensity factors k(a) and k(b) for sywnetric 
coll inear internal cracks in a strip. 


l<{a)/PQ/fc 


k(b)/Pp/a 
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APPENDIX A 


Exprssslons for tho functions K-j dnd ®nd the constents (see 
equation 17) (material type 1): 


, cosh(s,ax) 

= ?s^m,/( 5T ■^'"7 cosh(s,o.h) 

cosh(s2«x) 
+ 111^1111165) - nia cosh(s2<»hT 


+ m-jHigBg tanh(s^ah))3 > 


K2(x,a) 


cosh(s.^ctx) 
' 2 s 2 m, 3 P(a) ^"’7 cosh(s^ah) 


11 

12 

m, 


cosh(s2Ctx) 

2m^Q0g tanh(s2ah)) + mg cosh(S2ahI * 




P(a) = t3>^*i(s2ah) - m^mg tanh(s^ah) , 

m^ = l+^yx^l^l/^5 * "^2 ’ "^3 " 

m^ = S2+Vy/2 ’ "^5 " ''xy‘^'^1^1^^5 * "’e " ^5 ’ 

"*7 * * '"S ’’ W2‘^‘^2 ’ '"9 * 

"’lO ' ^2^2"”* ’ "’ll “ ’ "’12 ' ^2'h^h 


'"is ' ‘^r^2"’n^"’l2 * n>i4 = ('"5“'"6 "’i1^"’12^^^"’13 ’ 



*"l5 " '^3'"lo"V9 * ’ 

^17 ~ ■’^0&g(N^2^1 1*^i^li^g)/2Sinii 2 > 

in,8 = -ni7ro„8g(nt4+in2in,8/iii,2)/2m,3S2 , 

^19 “ nigHii ^ 3g(ni2'*’nigfn2/ni^ 2)/2s2nii 2 • 

APPENDIX B 

Expressions for the function K 2 and the constants r^ (see equa 
tion 20) (material type II); 

o 

K2(x,t,a) = y {[-rg Sin(w2ax)sinh(w^ax) 

+ rg cos((i)2ax)cosh(ujiax)3 • [r-jg sin[w2a(h-t)] • 

• (rg sin(w2ah)cosh(iD^ah) + r^g cosCwgoth) sinh(io-|ah)) 

- >'i8 

•(r^ sin(a)2ah)sinh((i>^ah) + r2 cos(w2ah) cosh(tUiah))3 


+ [rg sin(w2“^^s^*’^^^“l“^^ *^6 <^°s(u)2ax) cosh(w^cxx)3 '• 

• [-r^g sin[t»J2“(h-t)3 • (rgCOs(w2ah)sinh(u-|ah) 

jj 

- r^g sin(«2«h)cosh(w^ah)) + r-jg (-^ cos[w2ct(h-t)3 



- sin[t*J2a(h*t)]) • (-fg s1n(w2ati)sinh(u)iah) 


+ cosCwgCth) cosh(w^ah))]} , 

Q(a) = sinh(w^ah)cosh(w^ah) + rgQ sin{w 2 ah)cos((i) 2 (xh) . 

« w^+ iWg , $2 “ * ^^1^° ‘ 

c-j * * ^2 = * ^ 9^^ ^10 * ^2 ” ^ 9"^ ^10 * 

*"4 “ VS ^10‘‘'‘^2^9'^''^5 ’ *^5 ^ ’ ^6 “ ‘^2V^^8 ’ 

*’? ' ^xy‘*'^“l^9"“2^10^^^5 ’ ‘"S " V^yx ’ *'9 “ Wi^^-cOgPg-l , 

r^O ® ’ **11 " "^ 10 '*'‘^ 2 '^% ’ *^12 “ ’ 

*^13 ** “^9"^10'^12'^’"ll ’ **14 ' "^Wl2/'^n^/^’"l3 ’ 

**15 “ **4"*'3*’l2/'*ll ’ *^16 * **15V^*’l3'^l ’ **17 “ ■*'ir*'?2/*’ll 

**18 “ “l**l7/^*'l3 ’ **19 “ **1**10"*'2**9 * **20 ’*2**1 0’*’**! **9 ’ 

**21 ''^6*'l6**9"*'6*'l8*’l’*’*'5*'l0**16'‘‘*'2**5*'l8 ’ 

**22 = ^*'l**6“*’2*’5^'*18V'^l ’ 

**23 ' ■‘*'6*'l0*'l6"*'2*'6*'l8“*’5*’9*'l6"**l*'5**18 ’ 

*'24 = (V6'*’**l'*5^*'lgV‘^l ’ **25 "*'2r**24 ’ *’26 “ **2r*’24 

’*27 ” *'22"**23 ’ **28 ' **22'*'*’23 * **29 " ■(“i ''28‘‘’‘^2*'26^^^‘^l‘*’‘^2^ 




